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ABSTRACT 

Star-formation rates (SFRs) of galaxies are commonly calculated by convert- 
ing the measured Ha luminosities (Lna) into current SFRs. This conversion is 
based on a constant initial mass function (IMF) independent of the total SFR. 
As recently recognised the maximum stellar mass in a star cluster is limited by 
the embedded total cluster mass and, in addition, the maximum embedded star 
cluster mass is constrained by the current SFR. The combination of these two 
relations leads to an integrated galaxial initial stellar mass function (IGIMF, the 
IMF for the whole galaxy) which is steeper in the high mass regime than the con- 
stant canonical IMF, and is dependent on the SFR of the galaxy. Consequently, 
the Lhq-SFR relation becomes non-linear and flattens for low SFRs. Especially 
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for dwarf galaxies the SFRs can be underestimated by up to three orders of 
magnitude. We revise the existing hnear Lhq-SFR relations using our IGIMF 
notion. These are likely to lead to a revision of the cosmological star formation 
histories. We also demonstrate that in the case of the Sculptor dwarf irregular 
galaxies the IGIMF-formalism implies a linear dependence of the total SFR on 
the total galaxy gas mass. A constant gas depletion time scale of a few Gyrs 
results independently of the galaxy gas mass with a reduced scatter compared 
to the conventional results. Our findings are qualitatively independent of the 
explicit choice of the IGIMF details and challenges current star formation theory 
in dwarf galaxies. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
fundamental parameters — galaxies: irregular — stars: luminosity function, mass 
function — stars: formation 



1. Introduction 



Low-mass dwarf galaxies are very important objects to test the theoretical understand- 
ing of star formation, galactic evolution and cosmology. One key issue in this process is 
the accurate determination of current SFRs. Because massive stars have short life times 
the number of currently existing massive stars determines the actual massive star formation 
rate. The number of massive stars is commonly obtained from the integrated Ha luminosity 
of the target galaxy, after correcting the observations for extinction and N ll emission. If the 
fraction of massive star formation of the total star formation is known then the meassured 
integrated Ha luminosity of a galaxy can be converted into its current SFR. 

The standard method is to apply one initial mass function (IMF) to the entire population 
of newly formed stars in the whole galaxy. This standard method provides a linear relation 



betw e en the integrated Ha luminosity and the current SFR (lKennicuttlll983l : iKennicutt et al. 



1994: Kennicutt 



19981 ). Originally used for n ormal disk galaxies s uch linear relations have 



also been applied to dwarf irregular galaxies (jSkillman et al.ll2003l ). 



But even i f star formation in each young star cluste r follows the same canonical IMF 



( iKroupal 1200 ll . l2002l : iPfiamm-Altenburg &: Kroupal l2006l ) . the total stellar population of 

newly formed stars of all young star clusters together fo llows a distribution funct i on which is 

steep er than the canonical IMF in the high mass regime (IKroupa fc Weidneiil2003l : IWeidner fc Kroupa 
20051 ). This distribution function is called the integrated galaxial initial stellar mass func- 
tion (IGIMF) and deviates increasingly from the underlying canonical IMF with decreasing 
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galactic- wide star formation and total galaxy mass (IWeidner fc Kroupall2005l ). 



This effect has undoubtedly consequences for the cosmological evolution of the number 
of supernovae per low-mass star and the chemical enrichment of galaxies of different mass. 
As a first step the formulation of the IGI MF has been successfully used to n aturally explain 
the mass-met allicity relation of galaxies (IKoppen. Weidner fc Kroupal 120071 ) . 



Observational evidence supporting this IGIMF concept has been published by lSelman fc Melnick 



( I2OO5I ) who find a steeper IMF-slope of the high-mass regime in the field populati on of 30 Do- 
rad us tha n in th e cluster NGC 2070 as predicted by iKroupa fc Weidnerl (|2003). Note t hat 



the IScalol (119861 ) Milky- Way-disc IMF index a 2.7 flKroupa. Tout fc Gilmord Il993f ) is 
larger than the canon ical IMF index a fti 2.35 . This difference is a n atural consequence 
of the IGIMF notion jKroupa fc Weidnerl 120031 ) . Using the Scalo index IPiehl et al.l J2006h 
find a consistency of the supernova rate derived from the Galactic ^^Al gamma-ray flux 
and deduced from a surv ey from local O and B stars, extrapolated to t he whole Galaxy 
(iMcKee fc Williamsl ^9971 1 . Noteworthy is the finding bv llee et al] J2004h that LSB galax- 
ies have much steeper than Salpeter r nassive-star IMFs than M i lky- W ay-type galaxies, as 
expected from the IGIMF formalism. iHoversten fc Glazebrookl (120071 ) find evidence for a 
non-universal stellar IMF from the integrated properties of SDSS galaxies. Bright galaxies 
have a high-massive star slope of ~ 2.4, whereas fainter galaxies prefer steeper IMFs. The 
observ ed effect is very much in-line with the theoretical predictions of IWeidner &: Kroupa 
fl2005h . 



For a given global SFR the total number of massive stars based on an IGIMF is less 
than in the standard procedure. Thus, the expected IGIMF-Ha luminosity is less than the 
Ha luminosity calculated in the standard way, implying that SFRs of galaxies determined 
by the standard method are systematically underestimated. 

After a summary of the IGIMF basics we develop for different IGIMF-models a relation 
between the SFR and the produced Ha luminosity and provide easy-to-use fitting functions 
(Sec. [2]). In Sec. [3] we apply our Lh^-SFR relation to the dwarf irregular galaxies in the 
Sculptor cluster and discuss the changes of some of their parameters. Finally, we explore the 
detection limits for star formation (Sec. Hj). 



2. SFR-Lna-relation 

For a given SFR we assume that star formation takes place in star clusters distributed 
according to a mass function which is called the embedded cluster mass function (ECMF) 
of stellar masses (i.e. the function of the stellar-mass content in clusters before they expel 
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their residual gas), .Ccci(M:i)- These star clusters are formed within a time- span, 5t, which is 
called the "star-formation epoch" of the galaxy (jWeidner fc Kroupall2005l ). 

The total mass of all formed stars, Msfr, in all star clusters within the epoch 6t is then 
given by 

Msfr = SFR 6t = M,,i Ui{M,,i) dM,,i , (1) 

where Meci,min and Meci,max are, respectively, the least-massive and most-massive embedded 
cluster able to form. The total ensemble of all newly formed stars has a mass spectrum 
described by the IGIMF, C,igimf{^), constructed below. A star with the mass m produces 
a total number, Nion,st{m), of ionising photons during the epoch 6t. The total number of 
ionising photons emitted by all young stars within St is given by 

Nion,5t = / ^IGlMpim) Nion,Sti^) ' (2) 

where miow and mmax are, respectively, the minimum and maximum stellar mass. The 
calculation of A^ion,<5t will be described in Sec. 12.21 

Assuming that a fraction, /i, of all ionising photons leads to Ha emission in the sur- 
rounding gas then the resulting Ha luminosity is 

Lh„ = /i 3.0207 • 10-^2 erg A^io^/^t , (3) 

where the energy of one Ha-photon is 3.0207 ■ 10^^^ erg. Throughout this paper we assume 
/X = 1. It can be argued that the time-span 6t can be removed in this formulation by 
considering rates directly. But the total mass, Msfr, defined by Eq. [T] and depending 
explicitly on 6t, is an upper limit of the ECMF, as no cluster can be formed more massive 
than the available material. 



2.1. IGIMF 



LarsenI J2000l . l2002bl lal) and llarsen fc Richtleil J200oh found that the V-band luminosity 



of the b rightes t young cluster correlates with the global SFR. Based on the conclusion by 



LarsenI (l2002bl ) that the so -called super-clusters are just the young and massive upper end 



of a cluster mass function, IWeidner et al.l (120041 ) derived a relation between the maximum 
embedded cluster mass, Med max, in a galaxy and the current global SFR 



logio (Meci,max) = logiQ (/^ml) + 0.75 ■ log^Q (SFR) + 6.77 



(4) 
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where kuL is the mass-to-hght ratio, typically 0.0144 for young (< 6 Myr) stellar populations. 
They showed that this relation can be reproduced theoretically if the ECMF is a power-law 
with a slope of 2.35 and the entire star cluster population ranging from Mgci^min ~ 5 Mq to 
-^eci,max is bom withiu a time-span of 6t ^ 10 Myr, independently of the SFR. The concept of 
a relation between the SFR of a galaxy and the mass of the most mas sive star cluster has been 
devel oped further to derive the star formation history of galaxies (IMaschberger fc Kroupa 
2007h . 



Each cluster with the mass Med then forms stars between the mass limits mio 
'^max according to the canonical IMF, 



N 



and 



(5) 



This canonic al IMF, CM^^iim), is a two-part-pow er law and is here calculated using the 
algorithm by iPflamm-Altenburg fc Kroupal (120061 ). The normalisation constant of the IMF 



in an individual cluster, ^jvf^^j(m), with the mass Med an d the maximum stellar ma ss, m 



in this star cluster are determined by the two equations ( IWeidner fc Kroupal I2004J ) 



Mec\= I m ^Af,,i("^)dm , 

'miow 



(6) 
(7) 



The existence of an upper physical stellar of about 150 Mf7^ has been found in 

statistical examinations of several clusters ( Weidner fc Kroupal[200i : iFigeilbood : lOev fc Clarkl 
2OO5I : koenlbood : iMaiz Apellaniz et al.l boOfih. The canonical IMF. 



e 



km 



has the slopes 



as 



+1.30 
+2.35 
+2.35 



0.08 < tji/Mq < 0.50 
0.50 < m/Mo < 1.00 
1.00 < m/Mo < mma> 



(9) 



as used by IWeidner fc Kroupal (120051 ) who constructed different IGIMF-models. The third 
slope above 1 Mq has been introduced for varying the IMF-slope and constructing non- 
canonical models (as in Tab. [1]). 

The pair of the two integral expressions (Eq. 0, IT]) can not be solved analytically to 
obtain the relation between the star cluster mass , Mp^, and its most massive star, m me,.^, but 
the existence of a unique solution can be shown (jPflamm-Altenburg fc Kroupal 120061 ). Using 
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a simple bisection method this relation can be calculated numerically (Fig. [H crosses). The 
numerical solution can be very well fitted by the function 



y 



(10) 



with y = logio (mmax/M©), X = logio (M^ci/Mq), a = 2.56, b = 3.82, c = -0.38 and n = 9.17 
(Fig. m solid hue). 

This function behaves linearly for x b {y = (a/b) x + c) and is constant for x ^ 6 
(y=a+c). The transition occurs at x ~ 6 whereas a large n describes a rapid transition. 
Numerical simulations of star-formation in clusters (IBonnell et al.ll2004l ) also indicate that 
the mass of the most massive star scales with the system mass. 

The IG I MF i s then calculated by adding all stars in all clusters, as already noted by 



VanbeverenI fll982l . Il983[ ). 

^igimf("^) 



A^ecl,max('S'-F'-R) 



^Mcc\{m) ^ecl(Mecl) dM^d 



(11) 



M. 



Gcl,min 



2.2. Number of ionising photons 



The total number of ionising photons, A^ion,5t('^), emitted by a star of the mass m 
within the time 6t = 10 Myr, is calculated by integrating the emission rate of ionising 
photons, Nion{m,t), over 6t 



St 



N\on{m,t) dt. 



(12) 



Stars with life times shorter than 5t contributes all of their emitted ionising photons. 



In order to derive the emission rates of ionising photons, a grid of 4000 stars linearly 
spaced between 0.01 and 150 Mq is evolved from t o 20 Myr in time steps of 0.1 Myr. 
For stars below 50 Mq the stellar evolution package by lHurley. Pols fc ToutI (120001 ) is used. 
Above that limit u ntil 1 20 M©, the stellar parameters are interpolated from models by 
Meynet fc Maederl (120031 . without rotation ). For the most ma ssive stars above 120 JVIq 
formulae fitted to the massive star models by lSchaller et al.l ( 19921 ) are used in extrapolatior|^. 



The physical parameters provided by the models are then used to fin d appropriate stellar 
spectra for each star at ev ery time step from the BaSeL stellar library (ILejeune et al.l 119971 . 
19981 : IWestera et al.ll2002l ). From these spectra the part of the luminosity which can ionise 



^An extensive description of the fitting formulae can be found in lWeidner fc Kroupal (j2006l ) 
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logio(Mec|/Mo) 

Fig. 1. — The mmax-M^ci relation defined by Eqs. [6] and [7] with mmax* = 150 Mq and fitted 
by Eq. [10] (sohd curve). Only each 30th point of the numerical solution is marked by a cross. 
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Fig. 2. — The total number of ionising photons, N[on,st, emitted by a star with the mass m 
in solar masses within the time St — 10 Myr at the beginning of the hfe of the star. 
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hydrogen (all radiation with a wavelength A < 912 A), the ionising luminosity, Lion(^, 
is deduced for each star. Finally, the number of the ionising photons, Nion{m,t), can be 
calculated from 

iVion(m,t) = LionK^)10^o■^ (13) 
from Istahler fc Fallal J2005h . The resulting A^ion,a (Eq. ^ is shown in Fig. H 



In order to compare the results obtained with the here used set of stellar evolution mod- 
el^ with other models, the same procedure as d escribed above was appli ed to the Fadova94 
tracks for solar metallicity (IBressan et al.l Il993l ) which are preferred by iBruzual fc Chariot 
( 2003 ). As Bruzual &: CharlotI ( 20031 ) also use the BaSeL library of stellar spectra it is possi- 
ble to use our procedure to obtain A^ion but with the different stellar models. A^ion is plotted 
in Fig. [3] in dependence of stellar mass using the Padova94 tracks {dashed line) and using 
our models {solid line). Both compare very well. 



2.3. Models 

Although the IM F of young star clusters seems to be universal and widely independent 



of their environment (IKroupa 



2001 



20021). the universahty o f the ECMF slope is not an 



established result as summarised in IWeidner fc Kroupal (120051 ). To be consistent with our 
previo us work, the four IGIMF models here are the same as defined in IWeidner &: Kroupa 
koO^ . Both the IMF and the ECMF are multi-power laws, ^{m) oc m""' and ^eci(Meci) oc 



-ft 



eel 



All slopes and mass limits of all four IGIMF-scenarios are summarised in Tab. [H 



The Standard- IGIMF is based on a canonical IMF (Eq. [HD and an ECMF with a power 
P = 2.35 between 5 Mq and Meci,max(SFR). A slope of 2.35 provides the best agreement 
between the theoretical and the observed relation between the actual SFR and the mass of 



the heaviest young star cluster (IWeidner et al.ll2004l ). Furthermore, the theoretical formation 
time scale of star clusters has a constant value of 10 Myr independent of the cluster mass if 
a slope of 2.35 is chosen. 

In order to explore the dependence of the IGIMF on the ECMF we test three additional 
scenarios: 

The Minimal- 1 Scenario has a slope of /3i = 1 for small cluster masses (5 < M^ci/Mq < 
50) and /32 = 2 for higher masses (50 < Med/MQ < Med^max) and a canonical IMF (Eq. 



The Mimmal-2 Scenario has a truncated ECMF with p = 2 {50 < Mcd/M© < Med.n 



more detailed description of the models will be given in Weidner & Kroupa (2007, in preparation). 
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Fig. 3. — Logarithmic number of ionising photons (iVion,<ioMyr) versus stellar mass for the 
here used models {solid line) and the Padova94 tracks {dashed line). 
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Table 1. Slopes and mass limits of the IMF and ECMF for four IGIMF- models 



Parameter 


Standard 


Minimal- 1 


Minimal-2 


Maximal 




0.08 


0.08 


0.08 


0.08 




1.30 


1.30 


1.30 


1.30 




0.5 


0.5 


0.5 


0.5 


a2 


2.35 


2.35 


2.35 


2.35 


ma/Mo 


1.0 


1.0 


1.0 


1.0 


0:3 


2.35 


2.35 


2.35 


2.7 




5 


5 




5 


Pi 


2.35 


1.0 




2.35 


M2/M0 


50 


50 


50 


50 


02 


2.35 


2.0 


2.0 


2.35 



Note. — The IMF slopes aj refer to a power law, 
^{m) oc m~"', on a mass interval between mj and mi+i where 
^i+i < ^max*- The ECMF slopes Pi refer to a power- low, 
,^cci oc M~f', on a mass interval between M[ and Mj+i where 

Mi+i = Meel,^ax(SFR) (Eq. H. 
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i.e. no star clusters below 50 Mq are formed, and a canonical IMF (Eq. [9]). 

The Maximal Scenario consists of an IMF (Eq. ^ with a steeper high mass star slope 
of = 2.7 and an ECMF with (3 = 2.35 between 5 Mq and Meci,max(SFR). 

Both minimal scenarios produce a larger number of massive star clusters and therefore 
more massive stars than the Standard Scenario. The Maximal Scenario produces fewer 
massive stars than the Standard Scenario. Thus, for a given SFR both minimal scenarios 
produce a higher Ha luminosity and the Maximal Scenario a lower Ha luminosity than the 
Standard Scenario. 

In order to explore the change of the classical linear Lhq-SFR relation due to the effect 
of the IGIMF we compare the IGIMF models with two scenarios in which the IGIMF is 
invari ant, i.e. independent of the SFR. These invariant scenarios as well as the linear scenar- 



ios by iKennicutt et al.l (119941 ) provide theoretical Ha luminosities based on the assumption 
that the galactic-wide distribution function of newly formed stars corresponds to the IMF 
observed in young clusters. Thus, the slope of the IMF is independent of the total SFR and 
the IGIMF is identical to the IMF in young star clusters. 

In our first invariant model the IMF is chosen to have the canonical form defined by Eq. 
[9] and uses the same stellar properties calculated in Sec. I2.2l where mmax = ""^max* = 150 M0. 

Our second invariant model consists of a Salpeter IMF with a slope of 2.35 between 0.1 



and 100 Mq and the same stellar prope rties calculated in Sec 



the widely used linear model derived by IKennicutt et al.l (119941 ) 



221 This IMF is the basis of 



SFR(total) = ^f; nMoyr-^ • (14) 

1.26 • 10*^ erg s ^ 



2.4. Results 

The linearity of the widely used relation, Eq. [HI implies that if the total SFR is reduced 
by a certain factor then the produced Ha luminosity is reduced by the same factor, because 
the portion of ionising massive stars is constant if the IGIMF is identical to a universal 
IMF. But our concept of the IGIMF is a combination of two effects: firstly, with decreasing 
SFR the upper limit of the ECMF falls (Eq. Hj) and secondly the fraction of ionising massive 
stars is higher in heavy star clusters than in light star clusters due to the relation between 
the mass of a star cluster and its most massive star (Eq. [TOl) . Therefore, the galactic-wide 
fraction of massive stars is expected to decrease with decreasing total SFR. 

The resulting normalized Standard-IGIMF, i.e. the Standard-IGIMF divided by the 
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Fig. 4.— Normalised IGIMF in the Standard Scenario for different SFRs (10-^ 10-^ 
10~^, 10~^, 1, 100 yr~^), and the canonical IMF (upper-most hne). The arrow in- 
dicates the mass at which all IGIMFs start to deviate from the underlying canonical IMF 
("^max(Afeci,miii) = 1-25 M0). Notc that the IGIMFs are not exactly identical for m < 1.25 Mq 
due to the normalisation by the total mass of stars in each case. 
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total number of stars, is plotted in Fig. H] for different SFRs. As expected, all IGIMFs are 
steeper than the canonical IMF in the high-mass regime, and the mass of the most-massive 
star of the entire galaxy decreases with decreasing SFR. All IGIMFs start to deviate (Fig. HJ 
arrow) from each other above a mass threshold which is the maximum stellar mass of the 
smallest star cluster (mmax(Afeci,min = 5 Mq) = 1.25 M0). Stars below this threshold are 
formed in all star clusters. Note also, that the IGIMF has a lower upper mass limit than the 
canonical IMF. 

Because the functional form of the IGIMF is a function of the total SFR the produced 
Ha luminosity is expected to not depend linearly of the total SFR. The resulting relation 
between the SFR and the produced Ha luminosity is plotted in Fig. (5) This includes 
all four IGIMF-models and the two invariant models described abov e. Additionally, the 



widely used linear Lhq-SFR relation, Eq. [H] (IKennicutt et al.l Il994l . solid line), and the 



compl ete sample of linear rn odels based on different stellar evolution models and different 



IMFs (IKennicutt et al.lll994l . gray shaded area) are plotted, too. 



The IGIMF models are nearly linear above an Ha lumino sity of about 10 



40 



erg s 



which is the lower value of the normal disk galaxies explored by iKennicuttI (119831 ) , whereas 
the relations become much flatter below 10^° erg s^^. As expected the IGIMF-models lie 
above the linear relations. The order of magnitude of the deviation of the IGIMF-Lhq-SFR 
relation from the linear Lh^-SFR relation is independent of the explicit choice of the IGIMF 
model. 

For high Ha luminosities our invariant models are strongly linear like the Kennicutt 
relations. Below a SFR of about 10~^ M0 yr~^ they become shallower. In our invariant 
model the IGIMF is replaced by the canonical IMF. But the IMF is normalised such that 
the total mass contained in it is still given by 

SFR 5t= I m i{m) dm , (15) 



where mmax is the minimum of the physical upper mass limit, mmax* and the total mass 
content, SFR 5t (a star cannot be more massive than the total mass out of which it is 
formed). Note that ^(m) = ^uecxi.'^n) (Eq- ED for mmax = "^max*- If the mass available for 
star formation within the epoch 5t is less than the physical upper mass limit m^ax* then the 
upper limit of the IMF is given by the amount of available mass. Thus, at a SFR of 

SFR = ^ = = 10--«^ Mo (16) 

5t 10 Myr yr ^ ' 

our Lhq-SFR relations diverge from linearity. The simple physical boundary condition that 
stars more massive the SFR 6t cannot form leads to non-linear behaviour of the Lhq--SFR 
relation. 
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Fig. 5. — The calculated relation between the SFR and the resulting Ha-luminosity for four 
different IGIMF-models and two invariant models in which the IGIMF is assumed to be 
the canonical IMF and the Salpeter IMF . The solid line represents the widely used linear 



relation, Eq. [HI bv lKennicutt et al 



Lhq-SFR relations (IKennicutt et al. 



19941) . The gray shaded area marks the full set of linear 



1994 ) combining different stellar evolution models and 



different IMFs (for details see text). 
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For dwarf irregular galaxies having Ha luminosities between 10~^^ and 10~^^ erg the 
underestimation of the SFRs vary by a factor of at least 20 to 160 for the Minimal Scenario. 
In the Maximal Scenario the underestimation of SFRs vary between a factor of 800 and 5000. 
The SFRs of normal disk galaxies with an Ha luminosity of about 10^°-10^^ erg s~^ may be 
underestimated by a factor of 4 up to 50. 



2.5. Fitting functions 

For convenient use of the obtained Lhq-SFR relation we provide fitting functions for 
each IGIMF-model. The general form of the fitting function is a polynomial of fifth-order, 

y = a^x^ + a^x^ + a^x^ + a2X^ + aix + ao . (17) 

This function is used to fit the data-points of each IGIMF-model (Fig. [5]). The argument 
and the function value of the fitting-function are 

X = logio (^Ha [lO^^erg s"^] ) , 

/(x)=logio(SFR[Moyr-i]) . (18) 

The obtained values of the fitting parameters are listed in Tab. [2] for each model and the 
fitted models are plotted in Fig. [61 

The linear parts of our invariant models can be described by 

SFR = ^ M© yr-^ , 

1.89 • 10^1 erg s"! ° 

for Lh« > 2.8 • 10^^ erg s'^ (19) 
in the case of the canonical IMF and 

SFR = ,1;^ r Mo yr-^ , 

3.3 • 10*^ erg s ^ 

for LHa > 4.9 • 10^^ erg s"^ (20) 

in the case of the Salpeter- IMF. 



3. SFRs of dirr 



To study the effect of the IGIMF in real systems we now calc ulate the SFRs of the 
Sculptur dwarf irregular galaxies. The Ha luminosities are taken from lSkillman et al.l (120031 ) 
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Table 2. List of the fitting parameters of the Lna-SFR relation 



Model 


as 


04 


as 


a2 








L-Ra- 

(erg 


-range 

s-i) 




Standard 


-2.67e-05 


-9.30e-04 


-8.47e-03 


+3.21e-02 


+9.64e-01 


+4.38e-01 


1.7 


1025 


- 4.1 • 


10^3 


Maximal 


-2.94e-05 


-9.97e-04 


-9.13e-03 


+2.67e-02 


+9.63e-01 


+l.lle-00 


7.4 


1025 


- 1.2 • 


10^3 


Minimal- 1 


-2.32e-05 


-7.51e-04 


-5.70e-03 


+4.23e-02 


+8.85e-01 


-1.58e-01 


3.1 


1025 


- 2.2 • 


1044 


Minimal-2 


3.84e-05 


-2.43e-04 


-6.58e-03 


+3.50e-02 


+8.94e-01 


-1.91e-01 


8.6 


1032 


- 2.3 • 


1044 



Note. — Fit-parameters of the fitting-function (Eq. [T71) for the four different IGIMF-models: Standard, 
Maximal, Minimal-1 and Minimal-2 and the luminosity range of data points used for the fit. 




-10 -5 5 -10 -5 5 



logio(LHa/1o'^^ ergs-^) IoQiq (Lna / 10^^^ ergs"^) 

Fig. 6. — Fitted plot of the four IGIMF Lh^-SFR relations (crosses) using the fit function 
Eq. [17] and the parameters of Tab. [2] (solid lines). 
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and are already corrected for galactic extinction but not for N ll contamination as these 
dwarf galaxies have typically very low nitr ogen abundances. Th e SFRs are obtained using 



the fitting function Eq. f[T7|) . The SFRs in ISkillman et al.l (120031 ) are based on the classical 
relation Eq. (TH 

Observed Ha luminosities, linear and IGIMF-SFRs and H I masses, Mm , for eleven 
Sculptor dwarf irregular galaxies are listed in Tab. [31 The relation between the total galaxy 
gas mass, Mgas, and the SFR can be seen in Fig. [71 To account for Helium the total galaxy 
gas masses are calculated from 

Mgas = 1.32 Mhi , (21) 



following ISkillman et al.l fl2003r ) . 



The conventional data based on the Kennicutt relation (asterisks in Fig. [7|) can be 
divided roughly into two parts. Above a galaxy gas mass of about 2 ■ 10^ the SFRs are 
uniformly distributed around a value of about 7 ■ 10~^ M0 yr~^ with a large scatter ranging 
from 2.8 • 10"^ to 5.0 ■ 10"^ Mq yr'^ Below a galaxy gas mass of about 2 • 10^ Mq there is 
a steep decline of the SFR with decreasing galaxy gas mass. This relation is changed by the 
standard-IGIMF-effect (filled squares) in three ways: Firstly, all SFRs are increased. This 
increase is larger for low SFRs and therefore mainly for small galaxy gas masses. Secondly, 
the scatter of SFRs above a total galaxy gas mass of about 2 ■ 10^ Mq reduces by about one 
order of magnitude. Thirdly, the decrease of the SFR with decreasing total galaxy gas mass 
becomes less steep. 

A bivariate linear regression between the SFRs derived in the Standard-IGIMF model 
and the total galaxy gas mass (Fig. [3, solid line) leads to the relation 

'"^■"kSfT^I-OS- log. ^- 10.05 . (22) 

This means that the SFRs of the Sculptor dwarf irregular galaxies depends linearly on the 
total gas mass: 

SFR = 8.9- 10-^^ ^ -M^i^s . (23) 
yr ^ 6 

The corresponding gas depletion time scales. 



gas 



M,^/SFR , (24) 



are shown in Fig. [HI While the gas depletion times vary by a factor of 400 in the linear 
(Kennicutt) picture, the IGIMF leads to a shrinkage of the range down to only a factor of 
20. Additionally, the strong increase of the gas depletion time with lower galaxy gas mass 
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Fig. 7. — SFRs in Sculptor dwarf irregular galaxies in dependence of the total galaxy gas 
mass . The galaxy gas mass es and the SFRs calculated using the Kennicutt-relation are taken 
from I Skillman et al.l (l2003l . Tab. 3). These SFRs are compared with the ones obtained using 
the standard-IGIMF relation (SFRstd values in Tab. [2]) and with our invariant model where 
the IGIMF is replaced by a canonical IMF. The solid line marks the bivariate regression to 
the data calculated with the Standard-IGIMF described by Eq. 
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disappears and the gas depletion time becomes constant at a few Gyr. This is not surprising, 
as a strictly linear dependence of the SFR on the total gas mass, 

SFR = AM^^ , (25) 

directly implies a constant gas depletion time scale, 

T^^ = M^^/SFR = A-' . (26) 



4. Hct-invisible SF 



The determination of the SFR of a galaxy is based on the measurement of a total Ha 
luminosity of the galaxy and a reliable theoretical relation between the current SFR and the 
produced Ha radiation. The Ha emission is confined within distinct H ll-regions. These 
regions have to be identified manually after Ha imaging. The emission of these Ha regions 
have to exceed some detection limit to be distinguishable from background radiation and 
noise and to be detectable. 

Assuming that SF takes place in star clusters distributed according to an ECMF with 
star cluster masses down to a few solar masses then the least luminous observed H ll region 
is not the smallest star cluster but the cluster producing the lowest observable Ha radi- 
ation. Thus, the low est observed Ha flux indicates the detection limit. In the work by 
Skillman et al.l (l2003l ) this det ection limit is about 4 • 10~^^ erg cm~^ s~^ (H ll region ESO 
348-G9 No.2). In the work bv iMiller fc Hodgel Jl994h the lowest measured Ha flux in M81 



group dwarf galaxies is about 6.5 • 10~ erg cm~ s (IC 2574 MH 135). Assuming that 
each star cluster forms its own H ll region created by all ionising stars formed in this star 
cluster, a mass-luminosity relation (Fig. ^ between the embedded star cluster mass and the 
averaged produced Ha luminosity within the star formation epoch 6t can be constructed as 
described in Sec. 12.21 This can be well approximated by 



y 



{x + 34.55) (1 - 1/ exp(0.27x2 + 0.46x + 0.32)) 



(27) 



with y = logiQ(LHa/erg s ^) and x = \oR^ n(Mpc]/Mp)). For comparison t he Orion Nebula 
cluster with a total mass of abo ut 1800 Mrn (jHillenbrand fc Hartmannlll998l ) and an observed 



272000 Mq (jSelman et al.lll999l ) and an observed Ha luminosity of 1.5-10^° erg s ^ ( iKennicutt 
19841 ) are marked, lying well on the mass-Ha luminosity relation. 



Ha luminos ity of lO''^ erg s ^ (lKennicuttlll98^ and 30 Doradus with a total ma ss of about 



The Sculptor dwarf irregular galax y SC 24, for example , has only one observable H ll- 
region. Its flux is 8 ■ 10~^^ erg cm~^ s~^ ( jSkillman et al.ll2003l ) being just above the detection 
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Fig. 8. — Gas depletion times in Sculptor dwarf irregular galaxies in dependence of the total 
galaxy gas mass. The galaxy gas masses and the SF Rs calculated using the Lho-SFR relation 
by Kennicutt are taken from ISkillman et al.l (l2003l . Tab. 3). These SFRs are compared with 
the ones obtained using the Standard-IGIMF Lhq-SFR relation (Tab. [2]). 
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Table 3. SFRs of Sculptor dirr Galaxies 



Mm log;^o(-^Ha) SFRskiu SFRstd SFRmini SFRmin2 SFRmax 



Galaxy 


(10^ Mo) 


(erg s ^) 


(Mo 


yr 




(Mo 


yr 




(Mo 


yr 




(Mo 


yr 




(Mo 


yr 




ESQ 347-G17 


120 


38.9 


6.3 • 


10- 


-3 


4.1 • 


10- 


-2 


1.6 • 


10- 


-2 


1.4 • 


10- 


-2 


1.9 • 


10- 


-1 


ESQ 471-G06 


163 


38.4 


2.0 • 


10^ 


-3 


1.8 • 


10- 


-2 


7.9 • 


10- 


-3 


6.6 • 


10- 


-3 


8.0 ■ 


10- 


-2 


ESQ 348-G09 


84.3 


37.5 


2.8 • 


10- 


-4 


5.0 • 


10- 


-3 


2.6 • 


10- 


-3 


2.2 • 


10- 


-3 


2.1 • 


10- 


-2 


SC 18 


5.0 


36.5 


2.4 • 


10- 


-5 


1.5 • 


10- 


-3 


9.3 • 


10- 


-4 


8.4 • 


10" 


~4 


6.2 • 


10- 


-3 


NGC59 


16.7 


38.6 


3.1 • 


10- 


-3 


2.5 • 


10- 


-2 


1.0 • 


10- 


-2 


8.8 • 


10- 


-3 


1.1 • 


10- 


-1 


ESQ 473-G24 


63.8 


38.2 


1.3 • 


10- 


-3 


1.3 • 


10- 


-2 


6.0 • 


10- 


-3 


5.0 • 


10- 


-3 


5.8 • 


10- 


-2 


SC 24 


2.1 


35.4 


2.1 • 


10- 


-6 


5.8 • 


10- 


-4 


4.0 • 


10- 


-4 


4.0 • 


10- 


-4 


2.2 • 


10- 


-3 


DD0226 


33.9 


38.2 


1.3 • 


10- 


-3 


1.3 • 


10- 


-2 


6.0 • 


10- 


-3 


5.0 • 


10- 


~3 


5.8 ■ 


10- 


~2 


DD06 


15.4 


37.6 


2.9 • 


10- 


-4 


5.7- 


10- 


-3 


2.9 • 


10- 


-3 


2.4 • 


10- 


-3 


2.4 • 


10" 


-2 


NGC 625 


118 


39.8 


5.0 • 


10- 


-2 


2.2 • 


10- 


-1 


7.1 • 


10- 


-2 


6.3 • 


10- 


-2 


1.0 • 


10- 


-0 


ESQ 245-G05 


289 


39.1 


1.0 • 


10- 


-2 


5.9 • 


10- 


-2 


2.2 • 


10- 


-2 


1.9 • 


10- 


-2 


2.7- 


10- 


-1 



Note. — Integrated Ha luminosities and determined SFRs by ISkillman et al.l (120031 ) using the 
Kennicutt-relation (Eq. [T4l) . The SFRs of the dirrs for the IGIMF-models are calculated using our 
fitting functions. 
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limit. Given the distance of SC 24 of about 2.14 Mpc the Ha luminosity of the most luminous 
H ll-region in SC 24 is 2.5 ■ 10^^ erg s^^. Given our mass-luminosity relation (Fig. [9]) the 
corresponding star cluster stellar mass is about 280 Mq. Eq. IHthen determines the current 
SFR to be about 5.0 • 10~^ Mq yr~^ in good agreement with the Standard-IGIMF Scenario 
rather than the classical one (Tab. [3]). The classical SFR of 2.1 • 10~® M0 yr~^ implies 
that only 1 Mq of stellar material will have formed within 1 Myr which is by two orders 
of magnitude inconsistent with the expected cluster mass (280 M0). However, our SFR of 
5.8-10^^ means that 580 Mq must have formed within 1 Myr showing an internal consistency 
of our picture. 

An H ll-region with the luminosity Lhq- in a galaxy at a distance D produces an Ha 
flux of 

Using the mass-Ha luminosity relation (Eq. [27]) and a detection limit, jumit, a relation be- 
tween the star cluster mass and the maximum distance, -Dmax, within which the star cluster 
can be observed, can be constructed (Fig. [TO!) . 

If the Sculptor dwarf galaxy SC 24 were located more distant than 2.3 Mpc then SF 
in SC 24 would not have been detected. Therefore very moderate SF can be overseen. The 
effect of Ha-dark star formation on the cosmological SFR due to a possible large number of 
dwarf galaxies with low SFRs needs further study. 



5. Conclusion 

Using the integrated galaxial initial mass function (IGIMF) instead of an invariant IMF 
for calculating the pr oduced Ha-lumin o sity f or a given galaxy-wide SFR we revise the hnear 



Lhq-SFR relation by iKennicutt et al.l (Il994l ). The order of magnitude of the deviation of 
the IGIMF-Lna-SFR relation from the linear Lhq-SFR relation is independent of the explicit 
choice of the IGIMF model. 

Because the fraction of massive stars determined by the IGIMF is always less than given 
by the underlying invariant IMF of individual star clusters, SFRs are always underestimated 
when using linear Lhq-SFR relations. The linearity between the Ha luminosity and the SFR 
is broken, too. Therefore, the SFRs of normal disk galaxies with an Ha luminosity of about 
]^g4o_]^g42 g-1 j^g^y underestimated by a factor of 4 up to 50. For dwarf irregular 
galaxies having Ha luminosities between 10"'^^ and 10"'^^ erg s~^ the underestimation of 
the SFRs vary by a factor of at least 20 to 160 for the Minimal Scenario. In the Maximal 
Scenario the underestimation of SFRs of dwarf galaxies vary between a factor of 800 and 
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T 




1 2 3 4 5 6 

logio(Mec|/Mo) 

Fig. 9. — Mass-luminosity relation between the stellar mass of the cluster and the produced 
total Ha emission. The solid line shows the fit by Eq. [271 The observed position of the 
Orion Nebula cluster and 30 Doradus are indicated for comparison. Note the steep decline 
of for Meci < 1000 Mq is a result of the mmax(M5ci) relation (Fig. [T]). 
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Fig. 10. — Relation between the stellar mass of embedded star clusters and the maximum 
distance of the star cluster within which the star cluster can be detected with Ha observations 
in the case of three different detection limits: 10~^^, 10~^^ and 10~^^ erg s~^ cm~^. 
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5000. Evolutionary parameters of dwarf galaxies such as the gas depletion time scale are 
effected by the same order of magnitude. In fact, we find the IGIMF formulation to imply 
a linear dependence of the total SFR on the total galaxy gas mass and therefore the gas 
depletion time scale to be constant, a few Gyr, and independent of the total galaxy gas 
mass. 

If this finding is correct it imposes a problem for the theoretical understanding of galactic 
star formation explain ing the apparent low star formation efficiency of small galaxies (e.g. 



Kaufmann et al.l 120071 ). 



Observational evidence confirming the IGIMF scenario is becoming available (Hover- 
sten & Glazebrook, 2007) but further observational work is needed. Additionally, more 
consistency checks with different methods of deriving SFRs should be done in future. 

Further work is required to construct reliable conversions of IGIMF-fluxes into other 



pass bands. iKennicuttI (119981 ) mentioned that Ha measurements become unreliable for de- 
termining SFRs in regions with gas densities higher than 100 Mq pc~^ as present in star 
burst galaxies. An extension of the IGIMF-model to construct a relation between SFRs and 
FIR luminosities is therefore clearly necessary. 



We thank Ivo Saviane for stimulating discussions which encouraged us to write this 
paper. 



REFERENCES 

Bonnell I. A., Vine S. G., Bate M. R., 2004, MNRAS, 349, 735 
Bressan A., Fagotto F., Bertelh G., Chiosi C., 1993, A&AS, 100, 647 
Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 

Diehl R., Halloin H., Kretschmer K., Lichti G. G., Schonfelder V., Strong A. W., von Kienlin 
A., Wang W., Jean P., Knodlseder J., Roques J. -P., Weidenspointner G., Schanne S., 
Hartmann D. H., Winkler C., Wunderer C., 2006, Nature, 439, 45 

Figer D. F., 2005, Nature, 434, 192 

Hillenbrand L. A., Hartmann L. W., 1998, ApJ, 492, 540 
Hoversten E., Glazebrook K., 2007, ApJ, submitted (preprint) 



-27- 



Hurley J. R., Pols O. R., Tout C. A., 2000, MNRAS, 315, 543 
Kaufmann T., Wheeler C, Bullock J. S., 2007, ArXiv e-prints, 706 
Kennicutt Jr. R. C, 1983, ApJ, 272, 54 
Kennicutt Jr. R. C, 1984, ApJ, 287, 116 
Kennicutt Jr. R. C, 1998, ApJ, 498, 541 

Kennicutt Jr. R. C, Tamblyn R, Congdon C. E., 1994, ApJ, 435, 22 
Koen C, 2006, MNRAS, 365, 590 

Koppen J., Weidner C, Kroupa R, 2007, MNRAS, 375, 673 
Kroupa R, 2001, MNRAS, 322, 231 
Kroupa P., 2002, Science, 295, 82 

Kroupa P., Tout C. A., Gilmore G., 1993, MNRAS, 262, 545 
Kroupa P., Weidner C., 2003, ApJ, 598, 1076 
Larsen S. S., 2000, MNRAS, 319, 893 

Larsen S. S., 2002a, in Geisler D., Grebel E. K., Minniti D., eds, lAU Symposium Open, 
Massive and Globular Clusters - Part of the Same Family?, pp 421 — h 

Larsen S. S., 2002b, AJ, 124, 1393 

Larsen S. S., Richtler T., 2000, A&A, 354, 836 

Lee H.-c, Gibson B. K., Flynn C., Kawata D., Beasley M. A., 2004, MNRAS, 353, 113 
Lejeune T., Cuisinier F., Buser R., 1997, A&AS, 125, 229 

Lejeune T., Cuisinier F., Buser R., 1998, A&AS, 130, 65 

Mafz Apellaniz J., Walborn N. R., Morrell N. I., Niemela V. S., Nelan E. P., 2006, ArXiv 
Astrophysics e-prints 

Maschberger T., Kroupa P., 2007, MNRAS, 379, 34 

McKee C. F., Williams J. P., 1997, ApJ, 476, 144 

Meynet G., Maeder A., 2003, A&A, 404, 975 



-28- 

Miller B. W., Hodge P., 1994, ApJ, 427, 656 

Oey M. S., Clarke C. J., 2005, ApJ, 620, L43 

Pflamm-Altenburg J., Kroupa P., 2006, MNRAS, 373, 295 

Scalo J. M., 1986, Fundamentals of Cosmic Physics, 11, 1 

Schaller G., Schaerer D., Meynet G., Maeder A., 1992, A&AS, 96, 269 

Selman F., Melnick J., Bosch G., Terlevich R., 1999, A&A, 347, 532 

Selman F. J., Melnick J., 2005, A&A, 443, 851 

Skillman E. D., Cote S., Miller B. W., 2003, AJ, 125, 593 

Stahler S. W., Palla F., 2005, The Formation of Stars. The Formation of Stars, by Steven 
W. Stahler, Francesco PaUa, pp. 865. ISBN 3-527-40559-3. Wiley-VCH , January 
2005. 

Vanbeveren D., 1982, A&A, 115, 65 

Vanbeveren D., 1983, A&A, 124, 71 

Weidner C, Kroupa P., 2004, MNRAS, 348, 187 

Weidner C, Kroupa P., 2005, ApJ, 625, 754 

Weidner C, Kroupa P., 2006, MNRAS, 365, 1333 

Weidner C, Kroupa P., Larsen S. S., 2004, MNRAS, 350, 1503 

Westera P., Lejeune T., Buser R., Cuisinier F., Bruzual G., 2002, A&A, 381, 524 



This preprint was prepared with the A AS IM^jX macros v5.2. 



